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APPLIED ELECTROMAGNETICS 



THEORIE 



Direct and Related Far-Field Inverse 
Scattering Problems for Spherical 
Electromagnetic Waves in Chiral Media 

N.M. Berketis * 

Independent Researcher, Athens, Greece 

Abstract 

This paper studies the direct and inverse scattering 
problem when the incident electromagnetic field is a time 
harmonic point - generated wave in a chiral media and the 
scatterer is a perfectly conducting sphere. The exact 
Green's function and the electric far-field patterns of the 
scattering problem are constructed. For a small sphere, a 
closed-form approximation of the scattered wave field at the 
source of the incident spherical wave is obtained. Also 
treats the same inverse problem using far-field results via 
the leading order term in the low-frequency asymptotic ex- 
pansion of the scattering cross-section. 



Introduction 

In a homogeneous isotropic 
chiral media the electromag- 
netic fields are composed of 
left - circularly polarized 
( LCP ) and right - circularly 
polarized (RCP) components, 
which have different wave 
numbers and independent di- 
rections of propagation. 

The LCP and RCP components 
are assumed to be spherical 
Beltrami fields since in 
practice such wave fields are 
more readily realized. 

In this work, the author 
has studied the electromag- 
netic waves in chiral media 
produced by a point source in 
the vicinity of the scat- 



terer. In particular, in [1], 
[2], reciprocity, optical and 
general scattering theorems 
for stimulation of point- 
source asymmetric media have 
been demonstrated. This paper 
studies the inverse problem 
of far field [2] . Specifi- 
cally, we measure the scat- 
tering cross-section for a 
five-point source area. 

In the second Section, 
considering Bohren decomposi- 
tion into suitable Beltrami 
fields, we formulate the di- 
rect scattering problem of a 
spherical electromagnetic wa- 
ve by a perfectly conducting 
obstacle. This problem is 
well posed, the existence and 
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uniqueness has been proved in 
[1], [3]. 

In the third Section, af- 
ter expanding the incident 
field in terms of spherical 
wave functions, we obtain the 
exact solution of the scat- 
tering problem as well as an 
expansion for the electric 
far-field pattern [2]. 

Finally in the fourth Sec- 
tion, we consider either LCP 
or RCP incidence and we ob- 
tain an approximation of the 
scattering cross-section. For 
the far-field experiments, we 
measure the scattering cross- 
section for various point 
source locations [2] . 



vector r 0 with respect to an 
origin 0 in the vicinity of 
the scatterer. In order to 
define spherical electromag- 
netic fields E r , H, , we make 
use of the Bohren decomposi- 
tion into Beltrami fields 
Q t , r , and , as follows 

E r e = Qt,r s + Q»,r a 

1 ( 1 ) 

H, ~ (Ql-,r 0 " Q/?,r 0 ) 

Where q = (p/s) 1/2 is the in- 
trinsic impedance of the 
chiral medium. The Beltrami 
fields satisfy the equations 
[4], [5], 



Statement of the problem 

Our goal is to study the 
direct and inverse scattering 
problems when the incident 
electromagnetic field is a 
time harmonic point - gener- 
ated wave in a chiral medium 
and the scatterer is a per- 
fectly conducting sphere of 
radius a centered at the 
origin. The exterior space 
(r = |r| > a) is an infinite ho- 
mogeneous isotropic chiral 
medium with chirality measure 
(3, electric permittivity s 
and magnetic permeability p. 

We consider a time har- 
monic spherical electromag- 
netic wave due to a point 
source at P 0 with position 



Vx Ql.,r 0 ~ K<-Qf-,r 0 

^ X Q/?,r 0 = — V'rQi.,!-,) 



where y L andy R are wave num- 
bers given by, 



Yl l - kp ' Yr l + fcp 



with k = <y(sp) 1/2 , w being the 
angular frequency. The indi- 
ces L and R denote the LCP 
and RCP fields respectively. 
The spherical incident Bel- 
trami fields with suitable 
normalization have the form 
[ 1 ], [ 2 ], 
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Qi n r :(r|pJ = -(I + — VV + — Vxi) ( 

z y l y \ 



K (Yl |r ~r B |) 
K (y L r b) 



2nr g 1Yl ~ rQ ~ ~ 

£ [y L G fs (r, r 0 ) + VxG fs (r, r 0 ) ] • p L 



Q^(r|p R ) = -(i + - 5 -W- — Vxi) ( 
^ Yr Yr 



K (Yr |r -r B |) 
h 0 (y R r B ) 



2tt r e 1>v " e 

£ [y R G fs (r, r 0 ) - VxG fs (r, r 0 ) ] • p R 



(4) 



(5) 



where b 0 (x) = b 0 (x) = e lx /(ix) is 
the zeroth-order spherical 
Hankel function of first 

kind, I = xx + ijjijj + zz is the 
identity dyadic, r 0 =|r o | and 

G fs (r, r B ) is the free space 
dyadic Green function [4] . 
The constant unit vectors 
p L and p R satisfy the rela- 
tions 

V Pl= *V P>?= o 

r 0 xp L = ip, (6) 

r 0 xp R = -ip R 

We note that when r 0 -»oo, the 
incident electric field 

(T | p t , p w ) = 

( 7 ) 

= Qt"r" ( r I PJ + Qr™ (r I Pr) 



reduces to plane electric 
wave with direction of propa- 
gation -r 0 and polarizations 
P L , p R , since 

limQf-(r|pJ = e 1 ^ r p,= 

( 8 ) 

= Qf c (r;-f 0 ,p L ) 

limQ™(r |p R ) = e-^ r p R = 

r ^° ° (9) 

= Qf c (r; -r B , p R ) 

We consider E“ c is incident 

r 0 

upon a perfectly conducting 
sphere of radius a . Then, we 
want to calculate the scat- 
tered electric field E* c , 

which is the unique solution 
of the following exterior 
boundary value problem 
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( 10 ) 
( 11 ) 

The Silver-Muller radiation condition is modified as follows 
r x VxE* c (r) - Py 2 r xE* c (r) + E* c (r) = o(y ) , r ® (12) 



uniformly in all directions 
reS 2 , where S 2 is the unit 

sphere in M 3 , n is the out- 
ward normal unit vector on 

the scatterer and Y 2 =Y L Y R - 
The scattered electric field 
will be depended on the po- 
larizations p L , p R and will 
have the decomposition 

E r 0 C ( r I Pl , Pr) = Q.1% (r I Pz., Pr) 

. . (13) 

+ Q"r 6 (r I P„ Pr) 

where Q“ r< (r | p L , p R ) and 



The functions g L and g Rr 
are the LCP and RCP far-field 
patterns respectively [4], 
[ 6 ], 

If either a LCP or a RCP 
spherical electric wave 
Ef" c (r|pj, A = L, R , is inci- 
dent upon the scatterer, then 
the scattering cross-section, 
is given by [6], 



f lA |g,,r 0 (r|pj| 2 + 

s 2 Yl 

+ A |g R ,r 0 (r I pjf] ds{r) (16) 
Yr 



Q“r e (r I P„ Pr) are the corre- 
sponding scattered Beltrami 
fields which have the follow- 
ing behavior, when r oo 

^{r\p Ll p R ) = h.{ YL ry 

(14) 

•g L ,r 0 (r |p L ,p R ) + 0(^) 

Q"r e ( r IP^P R ) = h B (Y R r)- 

. . . (I 5 ) 

g R ,r 0 (r|P^PR) + 0(^) 



Exact Green's function 

We take spherical coordi- 
nates ( r,8,cp ) where 9 e [0, n] 
and (p e [0, 2tt ) , with the ori- 
gin at the center of the 
spherical scatterer, so that 
the point source is at 

r = r B , 0 = 0 . 



Thus, r 0 



= r 0 z , P,= 



72 »-**> 



and p R = —j= (x + iijj) , where x, ip 
v 2 



and z are unit vectors in 



FUNKTECHNIKPLUS # JOURNAL 



V2-10 



ISSUE 2 - YEAR 1 




DIRECT AND RELATED FAR -FIELD INVERSE SCATTERING PROBLEMS 



the x, l/j and z directions, 
respectively. Using spherical 
vector wave functions, [4], 
[5], [7], and taking into ac- 
count (4), and (5), we obtain 
[ 2 ], 

oo 

Qt,™ ( r I Pi.) = X { L oln(/l.r) + 

->=1 (17) 

+ iLW„(/ t r)} 



Q^(r|p R ) = Xe:{R-(/ R r)- 



-iRW(y R r)} 



where, for r < r B 



2n + 1 



" 2^/2h B (y A r B ) n(n + 1) 



H n ( y A r B ) 



Hn (W e ) = K ( y A r B ) - ih n ( y A r B ) 

with A = L, R . The h n is a 
spherical Hankel function of 
first kind of order n, 

h{x) = x -1 /7„(x) + h;(x) , 

and R^ n , with s = e, or 
o (even or odd) are the 
spherical functions [4], [5], 

^sL(yj) = M ii n (W r ) + N sln(W r ) # x 

( 20 ) 

RsiU^r) = H^,(y R r) - N^(y R r) 

where p = 1,3, the M^' n and 
N^ n are known spherical vec- 
tor function [7]. The scat- 
tered electric field that 
comes from a LCP incident 
field or a RCP incident field 
has a similar expansion to 
(17) or to (18) 



E“(r|pJ = Q“ (r|pJ + Q“,(r|pJ = 

00 oo 

= X B n a n { L o1„(W r ) + iL< e1n(/L r )} + X { R oln (K R 0 + (W) } 



Er 8 c (r I pj = Q“, (r | P R ) + Q ^ r0 (r | p R ) + 

00 00 

= X bX^Kyj) - + X 6X{R^(//) - ^£>*0} 



Using the boundary condition (11) on r = a, we obtain [2], 
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3n ( Y L a ) h n (y R a) + j„ (y L a) h n (y R a) 
h n (y L a) h n C y R a ) + h n {y L a) h n ( y R a ) 



in (V L a ) X - in. h n (/tg) 

h n (y t a) h n (y R a) + h n {y L a) h n (y R a) 



b L in (VVO X (ygg) - -7„ (WO h„ (y R a) 
h n (V'/.a) h n {y R a) + h n {y L a) h n {y R a) 



b R in (/gg) X Mi + j„ (y R a) /;„ (yya) (2g) 

" h n (y L a) h n (y R a) + h n (y L a) h n {y R a) 

Using the asymptotic forms [4], [7], 

L'si„ (y t r) - >(n + l)(-i) n b 0 (y,r) f s L ln (r ) ( 27 ) 

RsiU/rO ~ >(n + l)(-i) n b 0 (y R r) ff ln (r) ( 28 ) 

where let us introduce LCP Beltrami angular f s L ln (r ) , and RCP 
Beltrami angular f s R ln (r) [4], satisfy by relations, 

fsi n (r) = C sln (r) + iB sln (r) , f s R ln (r) = C sln (r) - iB sln (r) (29) 

we calculate the electric far-field patterns [2], 

sc .. . . . -Z (2n + l)(-i) n_1 W n (y,r 0 ) A A ... .... .... 

g a r (r p. ) = > - — , nV L 0/ a„ {fc„(r) - if 1n (r)} (30) 

Mxr e v iw Z. 2J2n(n + 1) b 0 (y L r o ) " x elnW olnW - r v ' 



sc ^ (2n + l)(-i) n 1 H n (y R r B ) A A ... ..„ .... 

g a r (r p.) = > - — . b„ {- f p 1 n (r) - if 1 n (r) } 

^ rR/ z. 2 l 2n{n + Jj hg{ y Rrg) nt el„U 0l „K J J 
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A far-field inverse problem 

So far, all of our formu- a 7 = - i 1 + ( y L a ) 3 + 0((y L af) 

las are exact. In the asymp- I 6 

totic results to follow, L . l + fik 5 7 

there are three parameters a 2 = - i — (Y l d) + 0((y L a)) 

y A a , with A = L, R and r = a/r 0 . 

We note that the geometrical y L a^Q (34) 

parameter r must satisfy 

0<r<l because the point r ifi-R/c') 3 

source is outside of the a i = - — ttttt (V L a Y + °((Y L a T) 

sphere. I A 1 + P K ) 

We assume that |y„a| « 1, as R i(l - (3/c) 4 

well; that is we make the so- 2 36(1 + (3/c) 3 L L 

called low-frequency assump- 
tion. From (23), (24), (26), y L a^O (35) 

(27) , we obtain [2] , 



, _ 1 + [3/c a 1 2n + 1 

" 2i£ 2 (2n + 1) 



R i (1 - (3/c) n + 2 , x2n + 1 



y L a 0 



bL ~ i (1 + p/c)- 2 
" 2n£ 2 (1 - (3/c) n + 1 1 



R i(l (3/c) y2n + l 
K 2C{2n + l) {YRa) 



blL= ~ 2(l-pff ^ a ) 3+ °(^ a ) 5 ) 

b 2 = - 1(1 + P/C) 3 (y R a) 5 + 0((y L a) 7 ) 
2 36(1 - (3/c) 3 * L 



y R a -> 0 



= - i(l T (3/c) (y R a) 3 + 0 ((y R a) 5 ) 
6 

(K R a) 5 +0((y R a) 7 ) 



y R a -> 0 



Vr 9 -> 0 ( 33 ) in order to calculate g*® 

where with an error of 0((y„a) 4 ) we 

only need the following [4] 

<■„= l-3-5 - (2n - 1) = (2n)\/(2 n n\) . 

In particular, 
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5(1 1 2^ T2[f - (?) “ ] > + 0((/,a) 4 ) 



where m A = 



-1, A = L 
1, A = R 



and if A = L,R at that case A C =R,L. 



Now the scattering cross-section, by LCP or RCP spherical 
Beltrami fields, from (16), is given by the relations [2], 

< c r e = \A |Cr s (r I P jf + 4 |g~, (r I P jf] cfs(r) = 



= (na 2 ){ (1 + |3/c) [j#r 4 



[f- r 4 + (y,a) 2 (M ^ + i 6 r 6 ) + (y L a) 4 (M + ^ r 4 )] 



(1 - (3/c) 4 3 



[- r + (y L a) (3r + 



5(1 - (3/c) 2 



(1 + (3/c) 2 4 16(1 + p/c) 2 



+ (Via) 4 (3 + 



4(3 + 9(1 + pff r4) ] } + °((^a) 6 ) , / t a -> 0 



0 *%, = f A fe ( ? I Pr )\ 2 + 4 fe (r I P R )| 2 ] rfs(f) = 



= (na 2 ){ (1 6 ^ /C) r 4 + (y R a) 2 (^ r 2 + 16r 6 ) + (y R a)\^- + ^ r 4 ) ] + 

(1 + (3/c) 4 3 4 , 2 ,_ 2 5(1 + (3/c) 2 

(1 - (3/c) l 4 ' K 16(1 -(3/c) 2 ' 



+ (/ R a) 4 (3 + 



4(3 + 9(1 -$k) 2 t4) ] } + °^ a ) 6 ) ' V rS 0 



In the special case r 0 ^- oo (r -+■ 0) , by the relations (46), 
(47) we obtain [2], 



. (1 + (3/c) 2 3(1 -(3/c) 4 

3 + (1 + p/c) 2 



} (Vi a ) 4 + 0((y L a) 6 ) , y L a -> 0 
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or 

cr = (tt a 2 ) (f«a) 4 + 0((/ fi a) 6 ) , y R a -> 0 (49) 

likewise in the case y L a ->■ 0 and y R a -+■ 0 , by the relations 
(46), (47) we obtain [2], 



< re = - f A {\ 3, /c) (na 2 ) (a/r 0 ) 4 



where A = L, R , with f,((3, /c) = 



(1 - ^(3/c) 2 ; 3(1 + ^(3/c) 4 
3 (1 - m A [ 3/c) 2 



Choose a Cartesian coordi- 
nate system Oxipz, and five 
point-source locations, name- 
ly (0,0,0), (1,0,0), (0,1,0) 

( 0 , 0 , 1 ) and ( 0 , 0 , 21 ), which 
are at (unknown) distances 
r 0 , r i i r 2 > r 3 and r 4 , respec- 
tively from the sphere's cen- 
ter. The parameter 1 is a 
chosen fixed length. For each 
location, measure the leading 
order term in the low- 
frequency expansion of the 
scattering cross-section. 

Thus, our five measure- 
ments are [2], 

fflj .= ^f 4 ((3, k)na\f-Y (52) 

J = 0, 1, 2, 3, 4 

Dimensionless quantities 
related to m j are 



v = _J_ = 2 ! 1 1 (H) 2 

Vj ^ if A (P,k)TT a ( a> (5J 

j = 0, 1, 2, 3, 4 

equivalently, we obtain [2], 



r] = | y'f,((3, /c)tt /j 
j = 0, 1, 2, 3, 4 



(54) 



There are six unknowns 
namely r 0 , r lt r 2 , r 3 , r 4 and a. 
Furthermore, r 0 , r 3 and r 4 are 
related using the cosine 
rule, r 2 + r 0 2 = 2r 3 2 + 2I 2 . So, we 
can find the six unknowns. 
The center of the spherical 
scatterer is obtained from 
the intersection of the four 
spheres centers at (0,0,0) 
(1,0,0), (0,1,0) and (0,0,1), 
with corresponding radius r 0 , 
r 1 , r 2 , r 3 respectively. 
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Abstract 

This paper investigates the influence of various parame- 
ters on the behavior of water droplets on polymeric surfaces 
under electric fields. An inclined plane test was carried 
out to understand the droplet behavior in strong electric 
field. Parameters such as, water droplet conductivity, 
droplet volume, polymeric surface roughness and droplet po- 
sitioning with respect to the electrodes were studied. The 
flashover voltage is affected by all aforementioned parame- 
ters. The droplet positioning is in some cases more vital 
than the droplet volume. 

Introduction face. Local arcing will even- 

Water droplets on a poly- tually bridge the dry bands 
meric surface may cause ar| d a complete flashover will 
corona under the influence of finally ensue. This mechanism 
an electric field and can is valid to a greater or 
cause deterioration to insu- lesser extent for both out- 
lation surface even in condi- door and indoor insulation, 
tions of low pollution level. although each of the afore- 
Water droplets on a polymeric mentioned categories have 
surface increase locally the their own particular charac- 
applied electric field. Local teristics, namely that indoor 
field intensifications lead insulation is stressed more 
to partial discharges (PD) and is subjected to a differ- 
and/or localized arcs, which ent type of environmental in- 
may render possible the dry fluences than outdoor insula- 
bands on the polymeric sur- tion [1], [2]. A combination 
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of water droplets and dust- 
like impurities on the sur- 
face of a polymeric surface 
may lead to a conducting con- 
tamination layer, which may 
cause a reduction of the 
flashover voltage. The design 
of high voltage insulators, 
they can be for indoor or 
outdoor use, one should take 
into account not only the 
pollution level, the insula- 
tor material and the appro- 
priate voltage level, but 
also the influence of water 
droplets on the flashover 
voltage. Previous work, car- 
ried out in this laboratory, 
tackled the behavior of water 
droplets for a wide range of 
water conductivities (1.7 
pS/cm - 10000 pS/cm) [3], [4]. 
In both publications, it was 
shown that, among the factors 
influencing the behavior of 
the water droplets, were the 
water conductivity, polymer 
surface roughness, droplet 
volume and droplet position- 
ing with respect to the elec- 
trodes . 

In the present work, a 
study of the aforementioned 
parameters on the water 
droplet behavior under the 
influence of a uniform elec- 
tric field in the range of 
1.7 pS/cm - 2000 pS/cm was 

carried out. All tests were 
performed with an inclined 
test arrangement, in order to 
simulate the behavior of wa- 
ter droplets on the surface 



of a real insulator. The an- 
gle used with respect to the 
horizontal was 10°. Such an 
angle was chosen because of 
its immediate relevance to 
industrial insulators. 

Force balance at the droplet/ 
polymer surface interface 

A modeling of a wet con- 
taminated surface was given 
in other publications and 
only a brief outline is pro- 
vided here [5] . Condensation 
of droplets on the surface of 
a high voltage insulator can 
come about from droplet 
germs. In Fig. 1, the forces 
exercised on the droplet are 
shown in case where no elec- 
trical field is applied. Such 
forces are the surface ten- 
sion of the liquid (t l ), the 
surface tension of the solid 
(t s ) and the interfacial ten- 
sion between liquid and solid 
(5 sl ) . When an electric field 
is applied, the droplet de- 
forms because of an addi- 
tional force. The tangential 
electric field on the surface 
of the insulator creates a 
force on the surface of the 
droplet which causes its de- 
formation. The deformation of 
the droplet affects the field 
distribution. Local field in- 
tensifications may result, 
which will cause micro-dis- 
charges between the droplets. 
This is the beginning of the 
chemical deterioration of the 
insulator surface. 
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Contact angle Water drop 




Fig. 1: Force balance at the interface solid/liquid at a 
water droplet on an insulating Surface (after [5]) 



Hydrophobicity may locally 
be lost. The voltage differ- 
ence across the droplet will 
be diminished and micro-dis- 
charges will follow. Solvable 
nitrates, which are the re- 
sult of the electrochemical 
deterioration, cause a higher 
conductivity of the water 
droplets. Dry zones may fol- 
low. It is important to bear 
in mind that not only the in- 
fluence of the applied elec- 
tric field on the shape of 
the droplet is of great sig- 
nificance, but also the in- 
fluence of the disintegrated 
droplet on the electric field 
distribution [5], [6]. 

Hydrophobic polymeric sur- 
faces are characterized by a 
low surface conductivity 
which in turn gives a low 
discharge activity and a 
higher flashover voltage. 
This holds also for polluted 
environments. Reduced hy- 



drophobicity implies a higher 
risk for flashover of the in- 
sulator. Hydrophilic materi- 
als, on the other hand, are 
very sensitive to polluted 
environments, and are charac- 
terized by a significant ac- 
tivity of local discharges 
[7] . It is well known the 
classification of the Swedish 
Transmission Research Insti- 
tute ( STRI ) regarding the hy- 
drophobic and hydrophilic 
surfaces. STRI classifies the 
various surfaces according to 
their hydrophobicity from 
class 1 (most hydrophobic, 
with only discrete droplets 
on the surface with contact 
angle larger than 80°) to 
class 7 (most hydrophilic, 
with continuous water films 
forming on the surface). How- 
ever, the truth is that no 
matter if the insulator has 
some sort of humidity and it 
is only slightly polluted or 
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it is heavily polluted, in 
both cases surface discharges 
play a most important role, 
and such discharges may start 
from water droplets. 

Experimental arrangement and 
preparation of the samples 

The aim of this paper is 
to study the behavior of wa- 
ter droplets under the influ- 
ence of an electric field. 
The voltage supplied was from 
a 20 kV transformer (in prac- 
tice the transformer may de- 
liver voltages up to 1.2 
times of its nominal voltage 
without loss of the accuracy 
of the measurement. Conse- 
quently, the applied voltages 
were accurate up to 24 kV) . 
The electrodes used were of 
copper. A top view as well as 



a cross section of an elec- 
trode is shown in Fig. 2. The 
electrodes were half cylin- 
drical in shape. Attention 
was paid to the smoothness of 
the electrode surfaces, so 
that no unnecessary field en- 
hancements could be noticed. 

The water droplets were 
positioned on the polymeric 
material surface with the aid 
of a special arrangement con- 
sisting of a metallic frame 
and three rules, one of which 
had two laser indicators. The 
water droplets were put on 
the surface with a syringe. 
Detailed information on the 
way the droplets were posi- 
tioned on the polymeric sur- 
face is given in [3] . The 
photograph of the inclined 
plane test is shown in Fig. 3. 





Fig. 2: Top view (above) and cross section (bottom) of the 
electrodes used (all dimensions in mm) 



FUNKTECHNIKPLUS # JOURNAL 



V2-22 



ISSUE 2 - YEAR 1 



EXPERIMENTAL RESULTS ON THE BEHAVIOR OF WATER DROPLETS 




Fig. 3: The inclined Plane test setup (side view) 



The polymeric materials 
used were PVC, rubber and 
silicone rubber. Surface 
roughness and resistivity of 
the material were measured. 
Surface roughness were mea- 
sured using perthometer (Type 
Perthometer M4P) . They gave a 
roughness of 0.25 pm for PVC, 
0.79 pm for silicone rubber 
and 1.10 pm for rubber. Re- 
sistivity of the material 
were performed with a Megger 
( BM25 type) and they gave a 
resistivity of 206 GQ for 
PVC, a resistivity of 3100 GQ 
silicone rubber and a resis- 
tivity of 2660 GO for rubber. 
The above values of surface 
roughness and surface resis- 
tivity were not isolated val- 
ues, but each of them was the 
mean of three measurements 
[8], [9], 

In the present work, by 



mixing known quantity of NaCl 
in distilled water forming 
solutions with conductivity 
in the range 1.7 pS/cm to 
2000 pS/cm, were used as 
droplet. The range of conduc- 
tivity were chosen based on 
the conductivity of natural 
rain and its values lie in 
the range 50 - 150 pS/cm, 
whereas the tests with porce- 
lain and glass insulators are 
performed with conductivities 
of 2500 pS/cm [10] . 

Experimental procedure 

The materials used were 
PVC, silicone rubber and rub- 
ber. Various droplet arrange- 
ments were studied. These ar- 
rangements are given in Fig. 
4. Each droplet had a volume 
of 0.2 ml. The electrodes 
were positioned at a distance 
of 4 cm from each other. 
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Fig. 4: Top view showing the droplet arrangements. Starting 
from top left, the arrangements were named as 
(1) arrangement with 1 droplet, 

(2a) arrangement with 2 droplets, 14-12-14, 

(2b) arrangement with 2 droplets, 12-16-12, 

(2c) arrangement with 2 droplets, 8-24-8, 

(3) arrangement with 3 droplets, 

(5) arrangement with 5 droplets and 

(9) arrangement with 9 droplets 

All dimensions given are in mm and they symbolize 
the distances of the droplets from the respective 
electrodes and the distances between them 



The parameters investi- 
gated were the water conduc- 
tivity, the roughness of the 
insulating surface, the posi- 
tioning of the droplets and 
their volume. The insulating 
surfaces were used as they 
were received from the manu- 
facturer without any further 
treatment. After putting the 



droplets on the polymeric sur- 
face, the voltage was slowly 
raised until flashover oc- 
curred. After that and after 
cleaning the surface and put- 
ting new droplets on it, the 
voltage was raised again up 
to the previous flashover 
value minus 1.2 kV, so that 
no new flashover would occur. 
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At this voltage the arrange- 
ment would stay for 1 min. If 
no flashover occurred, the 
voltage was raised by 0.4 kV 
and the procedure was re- 
peated until flashover oc- 
curred. The reason we left 
every time the voltage on for 
1 min, was in order to give 
necessary time interval for 
the droplet(s) to deform and 
for the partial discharge to 
initiate . 

It should be noted that it 
is observed a tendency for 
the droplets to slide, espe- 
cially for PVC because of its 
smooth surface. The droplet 
slide was minimal in the case 
of rubber, which was the 
rougher of the three materi- 
als used. An elongation of 
the droplets was observed, as 
the applied voltage was 
larger. A more evident oscil- 
lation of the droplet was ob- 
served with silicone rubber. 
The reason for that was be- 
cause the aforementioned ma- 
terial is more hydrophobic 
than the other two. Conse- 
quently, the droplet, for a 
defined droplet volume, has a 
smaller contact area with 
silicone rubber, and for this 
reason it oscillates more 
[8] . In some cases, such as 
with PVC with a droplet con- 
ductivity of 1.7 pS/cm and 
with the arrangement (1) of 
Fig. 4, ejection of minute 
charged droplets was observed 
just before flashover [11]. 



Experimental results 

At first, experiments were 
performed without any 
droplets between the elec- 
trodes. This was done in or- 
der to have reference values 
of the flashover voltage and 
also to understand influence 
of number of droplets between 
the electrodes that would re- 
sult in a reduction of the 
flashover voltage. The flash- 
over voltages without any 
droplets measured were 23 kV 
(± 0.5) for PVC, 25 kV (± 
0.5) for silicone rubber and 
24 kV (± 0.5) for rubber. The 
flash-over voltages of the 
three materials used were 
very similar. 

In Figs. 5 - 11 the varia- 
tion of flashover voltage 
with respect to the droplet 
conductivity for different 
droplet arrangements is 
shown . 

It is evident that sili- 
cone rubber presents a higher 
flashover voltage than the 
other two materials. It 
should be noted, however, 
that in the case of droplet 
arrangements (5) and (9) 
where rubber seems to be as 
good as silicone rubber. A 
possible explanation might be 
that in such a case, the 
droplets cover a significant 
part of the polymeric surface 
and hence they play an even 
more important role than the 
polymer itself. This in com- 
bination with the fact that 
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the rubber has a rougher sur- suit the lesser oscillation 
face compared to the other in the case of rubber, 
two materials, has as a re- 




Fig. 5: Flashover voltage for droplet arrangement (1) 




Fig. 6: Flashover voltage for droplet arrangement (2a) 
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Flashover voltage for droplet arrangement (2b) 




Fig. 8: Flashover voltage for droplet arrangement (2c) 
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Flashover voltage for droplet arrangement (3) 




Fig. 10: Flashover voltage for droplet arrangement (5) 
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Fig. 11: Flashover voltage for droplet arrangement (9) 



The better performance of 
silicone rubber is due to its 
hydrophobicity [5], [12]. The 

larger contact angle the 
droplets have minimum contact 
with the insulation material 
as in the case of silicone 
rubber. Figs. 12 - 14 show 

the influence of droplet vo- 
lume on the flashover volt- 
age. It is clear that the 
number of droplets affects 
the flashover voltage, i.e. 
the larger the number, the 
lesser the flashover voltage. 
An exception to that we have 
with the arrangement of 3 and 
5 droplets. It is evident 
that larger flashover volt- 
ages were observed with 5 
droplets than with 3 drop- 



lets. A possible explanation 
of that is that in the case 
of 3 droplets, the distance 
between electrode and droplet 
is smaller than in the case 
of 5 droplets. Consequently, 
one might say that there are 
occasions where the position- 
ing of the droplets with re- 
spect to the electrodes which 
plays a more vital role than 
the whole droplet volume. A 
further validation of the 
above consists of the compar- 
ison of the flashover volt- 
ages in the cases of 3 and 9 
droplets. It is observed that 
the flashover voltages for 
both these arrangements are 
not that different although 
the droplet volume triples. 
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An interesting case con- 
sists also of the droplets 
arrangements 2a, 2b and 2c. 
Higher values for flashover 
voltage were observed for ar- 
rangement 2a, then for ar- 
rangement 2b and the lower 
flashover voltage was ob- 
served for droplet arrange- 
ment 2c. This fact reinforces 
the above observations, name- 
ly that the positioning of 
the droplets play a crucial 
role, i.e. the closer the 
droplets to the electrodes, 
the lower the flashover volt- 
age. It is to be noted that 
similar observations were 
made also in [3], [4], where 

not an inclined arrangement 
was used but a horizontal 
one. What is presented in 
this paper is an approach of 
the behavior of water drop- 
lets on polymeric surfaces 
with an inclined electrode 
arrangement. The results were 
reproducible but not that 
many tests were carried out 
which would allow a statisti- 
cal study of the collected 
data. The main interest of 
this paper concentrates on 
the study of the behavior of 
the droplets. In the present 
context, no emphasis was gi- 
ven to the quantification of 
the studied parameters. 



A comment should be made 
on the results with the con- 
ductivity of 1.7 pS/cm: in 
such a case, with such a low 
conductivity, the water path 
behaves like a load, i.e. 
like a resistance connecting 
the two electrodes. It is for 
this reason that we observed 
a lowering of the voltage at 
the output of the power 
source. The flow of current 
through water of low conduc- 
tivity (i.e. of large resis- 
tance) means practically an 
increase of water temperature 
because of the power loss in 
the resistance of the water 
path. From the relation P = 
I 2 R, we can conclude that as 
the resistance is larger, as 
in our case the water path of 
very low conductivity, the 
power loss at this resistance 
is larger. Consequently, the 
temperature developed in such 
a resistance is enough for 
the boiling of the water. A 
quantity of water evaporates 
and the water path becomes 
narrower. Dry zones are 
formed, micro-discharges en- 
sue and finally the flashover 
follows. Such a phenomenon 
was observed in the inclined 
arrangement experiments, as 
they were observed before 
with non-inclined test ar- 
rangements [8] . 
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Conductivity (|jS/cm) 



Fig. 12: Flashover voltage for various conductivities, 
positioning and volume of the droplets. PVC used 




Conductivity (|jS/cm) 



Fig. 13: 
positioning 



Flashover voltage for various conductivities, 
and volume of the droplets. Silicone rubber used 
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Conductivity (pS/cm) 



Fig. 14: Flashover voltage for various conductivities, 

positioning and volume of the droplets. Rubber used 



Discussion and thoughts for 

further research 

In the present paper, some 
parameters influencing the 
droplet behavior on polymeric 
surfaces were investigated, 
such as water conductivity, 
droplet volume, polymeric 
surface roughness and droplet 
positioning. An increase of 
conductivity causes a de- 
crease of flashover voltage. 
This is a statement valid ir- 
respective of the polymer 
used. The surface roughness 
affects in a positive way the 
flashover voltage, when the 
number of droplets is large. 
The surface roughness func- 
tions as a hindrance to the 
movement of the droplets, and 
consequently renders their 
oscillation more difficult. 



An increase of droplet volume 
causes a decrease of flash- 
over voltage. This is in 
agreement with experimental 
observations published before 
with either ac or dc electric 
fields [13] . The position of 
the droplets with respect to 
the electrodes is of vital 
importance. With the droplets 
nearer the electrodes, the 
flashover voltage decreases. 
This is a phenomenon ob- 
served, albeit in different 
circumstances and conditions, 
also with enclosed cavities 
in solid dielectrics, where 
discharges become much more 
intense when one of the en- 
closing walls is an electrode 
[14]. 

The above show clearly 
that the polymeric material 
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plays a predominant role in 
determining the flashover 
voltage and the behavior of 
water droplets. Hydrophobic 
materials, such as silicone 
rubber, perform better than 
PVC or rubber. With this in 
mind, one should also note 
that most polymeric materials 
for outdoor applications pre- 
sent some sort of hydropho- 
bicity. However, the advan- 
tage of silicone rubber con- 
sists in the fact that it 
does not only have this prop- 
erty, it can also regenerate 
it [15]. 

The formation of water 
paths, between the droplets 
as well as between the drop- 
lets and the electrodes, gen- 
erally follow the direction 
of the applied electric 
field. The general activity 
in the form of discharges and 
droplet movement with rougher 
surfaces, sets in at higher 
voltages. In the case of just 
one droplet, with the appli- 
cation of the field, a defor- 
mation starts turning later 
to instability. Such behavior 
was observed with the in- 
clined arrangement as well as 
with previous horizontal ar- 
rangements [3], [4]. Also in 
the case of the inclined ar- 
rangement, the role of the 
'triple points' (i.e. the 
points where air, polymeric 
surface and droplet meet each 
other) is vital. The forces 
exercised on the droplets. 



because of the applied elec- 
tric field, are quite strong, 
and therefore, the 'triple 
points' move towards the 
electrodes. Experimental data 
published recently, validate 
what is reported here [16] . 
Such movement of 'triple 
points' causes the spread of 
the droplets. The spread of 
droplets is perhaps the most 
characteristic phenomenon ob- 
served with the inclined 
electrode arrangement. It is 
not, however, the only one 
observed. Droplet oscilla- 
tion, formation of water 
paths, collapsing of two 
droplets into a larger one, 
ejection of small charged 
droplets from a larger one, 
were also noted during the 
experiments. In this respect, 
the present work offers simi- 
lar conclusions with those in 
[3], [4], [8], [9], [16], [17], 
The importance of the triple 
points should be emphasized. 
In other works, it has been 
reported that partial dis- 
charge (PD) activity is some- 
times marginal, not easily 
detectable by a conventional 
electrical PD system but by a 
photomultiplier. Smaller wa- 
ter droplets offer higher 
flashover voltages. This can 
also be explained by consid- 
ering - in the case of smaller 
droplets - that the inter- 
molecular forces are in equi- 
librium with the surface ten- 
sion and, consequently, the 
electrical forces required to 



disturb the equilibrium are 
higher. On the other hand, 
for larger droplets, the in- 
ner forces may be lower and, 
if the surface tension stays 
the same, the equilibrium can 
more easily be disturbed. In 
such a case, the PD activity 
will start earlier and at 
lower voltages, and there- 
fore, the flashover voltage 
will be lower [18] . The fact 
that the contact angle in- 
creases with the decrease of 
water droplet volume, is 
something that cannot be un- 
derestimated [18] . Generally 
speaking, although the obser- 
vations of [18] refer to sil- 
icone rubber samples only, 
the relevant conclusions are 
not different from the ones 
presented in this paper. 

It is to be noted that the 
inclined electrode arrange- 
ment, used here, should not 
be compared by any means with 
the well known arrangement of 
the inclined plane test [19] . 
In the latter, a film of 
electrolyte is arranged to 
trickle down the back surface 
of a sheet and the samples 
are rated in terms of the 
voltage which causes a track 
to form in one hour [20] . In 
other words, the inclined 
plane test is a means of 
evaluating resistance to 
tracking and erosion of insu- 
lating materials for outdoor 
use, whereas the inclined 
electrode arrangement used in 



this work is a setup to study 
some parameters affecting the 
droplet movement on polymeric 
materials. The inclined plane 
test is an accelerating test 
[19] . The angle which was 
used in our experiments, i.e. 
the angle of 10°, was taken 
from real insulators. The 
purpose was to see the drop- 
let behavior under an elec- 
tric field in, as much as 
possible, real conditions. 
The present work confirms 
some general tendencies noted 
in [3], [4], 

It would be interesting to 
try experiments, in which the 
contact angle will be accu- 
rately measured for both 
smaller and larger water 
droplets w.r.t. time. The 
purpose would be to see 
whether the contact angle di- 
minishes with the passing of 
time and which diminution is 
more dramatic, that of the 
contact angle of smaller 
droplets or that of the con- 
tact angle of the larger 
droplets. Recent research in- 
dicated that the contact an- 
gle of droplets decreases 
with time, without, however, 
precising whether the rate of 
decrease is larger for smal- 
ler or for larger droplets 
[21] . Moreover, as noted be- 
fore, an interesting point 
can be to study different 
modes of droplet deformation 
in terms of the four parame- 
ters investigated in this pa- 
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per [18], [22], 

One last point should be 
raised: the research men- 
tioned above was carried out 
with conventional polymers. 
It is remarkable that, most 
of the aforementioned points, 
i.e. the importance of the 
positioning of the droplets, 
the significance of the drop- 
let volume and/or of the 
droplet number as well as the 
importance of the water con- 
ductivity, are points which 
are also important for non- 
conventional polymers, i.e 
for nanocomposite polymers. 
Work done recently in this 
laboratory emphasizes that 
parameters such as those in- 
vestigated here, are also 
significant for nanocomposite 
polymers. In fact, the nature 
of phenomena such as surface 
discharges or flashovers, must 
have a common underlying de- 
nominator for both conven- 
tional polymers and nanocom- 
posite polymers. Future re- 
search must also be directed 
in order to find the common 
grounds for surface discharge 
phenomena in both conven- 
tional and nanocomposite poly- 
mers [23] , [24] . 

Conclusion 

Water droplet conductiv- 
ity, polymer surface rough- 
ness, droplet volume and the 
positioning of droplets with 
respect to the electrodes 
constitute important parame- 



ters affecting the behavior 
of droplets under the influ- 
ence of an electric field 
with an inclined plane elec- 
trode arrangement. Increased 
conductivity, smoother poly- 
mer surfaces and increased 
droplet volume cause a reduc- 
tion of the flashover volt- 
age. The droplet positioning 
with respect to the elec- 
trodes plays a vital role in 
reducing the flashover volt- 
age and, on occasions, is 
more important than the drop- 
let volume. 

Appendix - Elementary 

modeling 

As mentioned above, the 
behavior of a droplet was 
modeled in [25], where the 
electric field E N developed 
at one edge of the droplet is 
given by 

E N = U h/[a(h - a)] (1) 

where, U is the applied volt- 
age, a is the droplet radius 
and h is the distance of the 
center of the droplet from 
one of the electrodes. The 
electric field E„ on the 
opposite edge of the droplet 
is given by 

E M = Uh(L-h)/[a(L-h-a)] (2) 

where, L is the distance be- 
tween the electrodes and the 
other symbols as in Eq . (1). 
As a droplet is positioned in 
the middle of the electrodes, 
the ratio E N /E M is expected to 



be unity, and this was what 
was exactly obtained with the 
above equations. Considering 
Eq . (1), as h^a, E N tends to 

infinity. This is what we 
observed in the context of 
this work, namely that the 
closer a droplet is in one of 
the electrodes, the larger 
the electric field is get- 
ting, and consequently the 
more deleterious the conse- 
quences are. The above simple 
modeling is due to [25] . It 
was elaborated in [26], 
albeit with a non-inclined 
test arrangement. It must be 
mentioned that phenomena of 
very similar nature were 
observed also with the in- 
clined test arrangement. 



Remark 

Main aspects of this paper 
were published in M.G. 
Danikas, P. Ramnalis and R. 
Sarathi, "A study of the be- 
havior of water droplets on 
polymeric surfaces under the 
influence of electric fields 
in an inclined test arrange- 
ment", Journal of Electrical 
Engineering, Vol. 60, No. 2, 
2009, pp. 94-99. The present 
version, however, contains 
additional comments on some 
aspects of the investigated 
polymeric materials as well 
as on a possible relation of 
the observed phenomena in 
nanocomposite polymers. 
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Abstract 



In order to demonstrate the usefulness of the only one 
existing method for systematic error estimation in VNA (Vec- 
tor Network Analyzer) measurements by using complex DERs 
(Differential Error Regions), we compare one-port VNA mea- 
surements after the two well-known calibration techniques: 
the quick reflection response, that uses only a single S 
(Short circuit) standard, and the time-consuming full one- 
port, that uses a triple of SLO standards (Short circuit, 
matching Load, Open circuit). For both calibration tech- 
niques, the comparison concerns: (a) a 3D geometric represen- 
tation of the difference between VNA readings and measure- 
ments, and (b) a number of presentation figures for the DERs 
and their polar DEIs (Differential Error Intervals) of the 
reflection coefficient, as well as, the DERs and their rec- 
tangular DEIs of the corresponding input impedance. In this 
paper, we present the application of this method to an AUT 
(Antenna Under Test) selected to highlight the existence of 
practical cases in which the time consuming calibration tech- 
nique results a systematic error estimation stripe including 
almost all of that of quick calibration. 



Introduction 

The systematic error in a 
full one-port calibrated VNA 
measurement p of a given one- 
port DUT (Device Under Test) 
is already estimated by its 
DER [1] - [2] : 

p = (m - D)/[M(m - D) + R] (1) 



dp = [- RdD - (m - D) 2 dM - 
- (m - D)dR + Rdm] / 

/ [M(m - D) + R] 2 (2) 

where m is the VNA complex 
reading and D, M and R are 
the complex system errors of 
Fig. 1. 

The relations holding be- 
tween this complex reflection 
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coefficient p and its respec- 
tive impedance Z, as well as, 
between their DERs are [1] - [2] : 

Z = Z 0 (l + p)/(l-p) (3) 

dZ = 2Z 0 dp/(l - p) 2 (4) 

In this paper, we express 
the DERs for systematic error 
estimation in VNA measure- 
ments calibrated by the much 
simpler and quicker reflec- 
tion response technique, in 
order to be in place to make 
some practical decisions from 
the different calibration tech- 
niques comparison. 



1 




R 



Fig. 1: Full one-port error 
model 

Response Calibration 

The reflection response ca- 
libration technique can be 
accomplished with the mea- 
surement of only one standard 
load, instead of three in 
full one-port, usually of a S 
short circuit [3] -[4]. This 
means that the flow graph of 
Fig. 1 is simplified a lot, 
since the two surrounded by 
dashed boxes system error 
branches of directivity D and 
source match M do not exist, 
equivalently D = 0 and M = 0 
and (1) results to: 



R = m/ps = s/S (5) 

where s is the VNA complex 
reading of the S short cir- 
cuit standard with a nominal 
value of S = - 1, m is the com- 
plex reading of a given DUT 
and ps is its complex reflec- 
tion coefficient as it is 
measured after this response 
calibration : 

p s = (m/s)S (6) 

which, from (2), has the dif- 
ferential error: 

dps = (S/s)dm - (Sm/s 2 )ds + 

+ (m/s)dS (7) 

The corresponding total DER 
is then the sum of L = 3 par- 
allelograms. Therefore, this 
DER contour is a polygonal 
line with 4L = 12 line seg- 
ments and vertices at most, 
in contrast with the DER of 
the measurement after a SLO 
full one-port calibration, 
which is a piecewise curve com- 
posed of 4(L - 1) = 24 line 
segments, 4(L - 1) = 24 circu- 
lar arcs and 8(L-1) = 48 ver- 
tices, at most [1] - [2] . 

Application Results 

By following the error es- 
timation process, we already 
detailed in [l]-[2], we take 
as dS the considered manufac- 
turers' standard S uncer- 
tainty data: 

-0.01 < d|S| < 0, -2° < dzS < +2° 



FUNKTECHNIKPLUS # JOURNAL 



V2-42 



ISSUE 2 - YEAR 1 



VNA UNCERTAINTY PART 3: SHORT ONE-PORT CALIBRATION - COMPARISON 



and as dm and ds the VNA in- 
accuracy of ±1 digit in LSD 
of their corresponding read- 
ings, for either the ampli- 
tude in decibels or the phase 
in degrees. Moreover, the one- 
port DUT that was considered 
is the same typical UHF ground- 
plane antenna (that is: AUT) 
mentioned in [1] . 

The difference between the 
3 nominal values (-1, 0, 1) 
of the 3 full one-port cali- 
bration standards (S, L, 0), 
respectively, and their 3 cor- 
responding VNA readings (s, 

1, o), can be estimated by 
the extent of the surfaces 
shown in the triptych of Fig. 

2, where the vertical axis 
segment represents the range 
of the distinct stepped fre- 
quencies. Each surface is for- 
med by parallel to horizontal 
plane lines. Each such line 
expresses the complex differ- 
ence between the standard no- 
minal value and its corre- 
sponding VNA reading, in each 
stepped frequency. 

In the triptych of Fig. 3, 
and from left to right we 
have the difference surfaces 
made by distance lines be- 
tween : 

(a) the measured reflection 
coefficient p after a full 
SLO one-port calibration 
(black solid points) and the 
corresponding VNA readings m 
for the AUT measurement (col- 
ored magenta points). 



(b) the measured reflection 
coefficient ps after S re- 
sponse calibration (black 
ring points) and the corre- 
sponding VNA readings m for 
the AUT measurement (colored 
magenta points), and 

(c) the two measurements (p, 
Ps) . 

All the involved, previ- 
ously shown, quantities are 
projected on the horizontal 
complex plane of Fig. 4. The 
magenta colored spiral repre- 
sents m, while, the black 
curves the reflection coeffi- 
cient: solid points, for p, 

and ring points for ps. All 
of 1 VNA readings are close 
enough to complex origin (co- 
lored green points). It is 
rather difficult to distin- 
guish the two curves for s 
and o VNA readings, which are 
close enough to the unit cir- 
cle circumference (colored 
red solid points and colored 
blue ring point, respec- 
tively) . 

The p-DERs and ps-DERs, 
for all 4 MHz stepped fre- 
quencies covering the range 
of [600, 1000] MHz, are over- 
lapped on the complex plane 
of Fig. 5, forming a light 
and a dark gray stripes, re- 
spectively. From each stripe 
we selected 11 DERs out of 
101, drawn with dark gray and 
white colors respectively, to 
illustrate their outline de- 
pendence on frequency. 
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Fig. 4: VNA s, 1, o, m readings and p, ps measurements. 



Moreover, we selected to 
magnify a part of this figure 
in the sub-range of [892, 
1000] MHz, to further illus- 
trate the DER outlines and 
their overlapping in Fig. 6, 
where the clearly shown rip- 
ple of the simple response 
calibration stripe over the 
relatively smooth full one- 
port calibration stripe re- 
veals the superiority of the 
latter in the production of 
more accurate measurements. 

The comparison between the 
AUT measurements based on 



these two calibration tech- 
niques is extended to the 
comparison against the fre- 
quency : 

(a) of the computed polar 
DEIs of the reflection coef- 
ficient magnitude and argu- 
ment stripes in Fig. 7, 

(b) of the rectangular DEIs 
for the corresponding R input 
resistance and X input reac- 
tance stripes, in Fig. 8 and 

(c) of the Z-DERs, and Zs- 
DERs stripes in Fig. 9. 
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Fig. 5: Complex p-DERs and ps-DERs in [600, 1000] MHz 




Fig. 6: Complex p-DERs and ps-DERs in [892, 1000] MHz 
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Fig. 8: Rectangular DEIs of input impedance 



SUNDAY 15 DECEMBER 2013 



V2-47 



FUNKTECHNIKPLUS # JOURNAL 



Input Reactance X, Ohm 3 Reflection Coefficient Phase, [°] 



N.I. YANNOPOULOU, P.E. ZIMOURTOPOULOS 




Fig. 9: Complex Z-DERs and Zs-DERs in [600, 1000] MHz 



Conclusion 

From all that, it must be 
clear now that in this inten- 
tionally selected for presen- 
tation particular AUT case 
there was no advantage at all 



in selection of full one-port 
calibration over the reflec- 
tion response one, due to their 
remarkable in all aspects co- 
incidence. Of course this is 
just another one conclusion 
a-posteriori . 
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A faded synthesis of an anthemion rooted in a meandros 

The thirteen-leaf is a symbol for a life tree leaf. 

"Herakles and Kerberos", ca. 530-500 BC, 
by Paseas, the Kerberos Painter, 
Museum of Fine Arts, Boston. 

www. mfa . org/collect ion s/object /plate -153852 

The simple meandros is a symbol for eternal immortality. 

"Warrior with a phiale", ca. 480-460 BC, 

by Berliner Maler, 

Museo Archeologico Regionale "Antonio Salinas" di Palermo, 
commons .wikimedia . org/wiki/File :Warrior_MAR_Palermo_NI2134. j pg 
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